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ABSTRACT 


This report provides an overview of the Hyper-X research vehicle Monte Carlo analysis 
conducted with the six-degree-of-freedom simulation. The methodology and model uncertainties 
used for the Monte Carlo analysis are presented as permitted. In addition, the process used to 
select hardware validation test cases from the Monte Carlo data is described. The preflight Monte 
Carlo analysis indicated that the X-43A control system was robust to the preflight uncertainties and 
provided the Hyper-X project an important indication that the vehicle would likely be successful 
in accomplishing the mission objectives. The X-43A in-flight performance is compared to the 
preflight Monte Carlo predictions and shown to exceed the Monte Carlo bounds in several instances. 
Possible modeling shortfalls are presented that may account for these discrepancies. The flight 
control laws and guidance algorithms were robust enough as a result of the preflight Monte Carlo 
analysis that the unexpected in-flight performance did not have undue consequences. Modeling 
and Monte Carlo analysis lessons learned are presented. 

NOMENCLATURE 


alpha 

angle of attack 

AIL 

aircraft-in-the-loop 

AMW 

all-moving wing 

Avg 

average 

CA 

axial force coefficient 

CA bias 

axial force bias 

CFD 

computational fluid dynamics 

CG 

center of gravity 

Cl 

roll moment coefficient 

Cl bias 

roll moment bias 

Clp 

roll moment caused by roll rate 

Cm 

pitch moment coefficient 

Cm bias 

pitch moment bias 

Cmo 

pitching moment at zero degree angle of attack 

Cmq 

pitch moment caused by pitch rate 

Cn 

yaw moment coefficient 

Cn bias 

yaw moment coefficient bias 

Cnp 

yaw moment caused by roll rate 

Cnr 

roll moment caused by yaw rate 

CP 

pressure coefficient 

CY 

side force coefficient 


CY bias 

side force bias 

CZ 

normal force coefficient 

CZ bias 

normal force bias 

DFRC 

Dryden Flight Research Center 

FADS 

flush airdata sensing 

FMU 

flight management unit 

G 

gravitational acceleration 

GNC 

guidance, navigation, and control 

GRAM 

Global Reference Atmospheric Model 

h 

altitude 

HIL 

hardware-in-the-loop 

HXLV 

Hyper-X Launch Vehicle 

HXRV 

Hyper-X Research Vehicle 

INS 

inertial navigation system 

kft 

thousands of feet 

MCAT 

Monte Carlo analysis tool 

PID 

parameter identification 

PPT 

pressure port transducer 

SepSim 

detailed simulation of separation dynamics 

SIM 

simulation 

UA 

unavailable 

6-DOF 

six degree-of-freedom 

a 

drag dispersion 


INTRODUCTION 

Before the initial flight-testing of a new aircraft, there is always some uncertainty as to how 
the vehicle will actually perform. These uncertainties increase when the vehicle is a unique aircraft 
flying in a seldom-explored flight regime. One such unique vehicle is the X-43A, which was 
the first unmanned vehicle to fly at hypersonic speeds with a scramjet engine. A Monte Carlo 
analysis was conducted to predict and bound the performance of the X-43A, hereafter referred to 
as the Hyper-X research vehicle (HXRV). Monte Carlo analysis is a useful and commonly used 
technique to assess system performance to predicted parameter variations by means of a stochastic 
process. The analysis consists of independent and random variations of all applicable model 
inputs and parameters within the estimated uncertainty bounds. Random combinations of input 
uncertainties can produce different, and potentially worse, results than deterministically varying 
a single uncertainty or set of uncertainties. Monte Carlo analysis consists of making multiple 
simulation runs, with each run containing a unique set of parameter variations, and analyzing the 
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results. A close examination of these simulation runs provides insight into the system performance 
and characteristics. Monte Carlo analysis forces the flight control laws and guidance algorithms 
to be robust to a large range of uncertainties. Stress testing past the modeled uncertainty range can 
provide an assessment of the vehicle robustness to the unknown unknowns present in every flight 
research program. This type of analysis is particularly useful for unmanned vehicles where there is 
no pilot in the loop acting as an adaptive controller. 

The purpose of this report is to provide an overview of the HXRV Monte Carlo analysis tools, 
present the uncertainties used as part of the Monte Carlo analysis (as permitted), and discuss the 
results of the Monte Carlo analysis. In addition, limited comparisons between simulation and flight 
data are presented with an emphasis placed on highlighting modeling shortfalls. The selection 
process for the subset of runs used during hardware validation testing is described, and lessons 
learned and observations concerning the Monte Carlo analysis are presented. 

PROJECT DESCRIPTION 

The goal of the NASA Hyper-X program was to demonstrate an advanced, airframe- 
integrated, airbreathing, hypersonic propulsion system in flight, thus validating the supporting 
tools and technologies (ref. 1). The program sought to design, build, and flight-test a series of three 
small, autonomous, scramjet-powered Hyper-X research vehicles (HXRV), designated X-43A, 
at Mach 7 and 10. The vehicles were dropped from the NASA Dryden Flight Research Center 
(DFRC) B-52B airplane (The Boeing Company, Chicago, Illinois) over the Pacific Ocean and 
rocket-boosted to its test point. Prior to flight test, the Hyper-X program conducted extensive 
ground testing, simulation, and analysis. In support of these activities, a high-fidelity six degree- 
of-freedom (6-DOF) nonlinear simulation of the research vehicle was developed at DFRC. The 
research vehicle simulation was meant to accurately model the vehicle dynamics after separation 
from the Hyper-X launch vehicle (HXLV), including the engine test and controlled flight down to 
the water. Many of the simulation models included uncertainties or dispersions that were used in 
the Monte Carlo analysis. 

The Hyper-X Research Vehicle (HXRV) Description 

As shown in figure 1, the HXRV was an unmanned autonomous vehicle that measured 
approximately 12 ft long and 5 ft wide, and weighed approximately 3000 lb. The scramjet engine 
was attached to the underside of the HXRV and contained a cowl door on the leading edge of the 
engine that was used to control airflow through the engine. The HXRV had four control surface 
effectors: a left and right all-moving wing (AMW) and twin rudders. The HXRV was fitted to 
the front end of the Hyper-X launch vehicle (HXLV), a modified Pegasus® (Orbital Sciences 
Corporation, Chandler, Arizona) rocket used to boost the HXRV to the desired test condition. 
The combination of the HXRV and HXLV was known as the X-43A stack, and was carried to the 
launch point under the wing of the NASA DFRC NB-52B airplane, ship number 008. 
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Figure 1. The X-43A 3-View. 


Flight Test of the Hyper-X Research Vehicle 

Three HXRVs were built for the Hyper-X program. Two of the vehicles were intended 
for missions at Mach 7, and one vehicle was intended for a mission at Mach 10. The HXRVs 
collectively received the designation of X-43A. All three vehicles have the same outer mold line; 
the primary difference among the vehicles is the internal engine flowpaths. Flight 1 was intended 
to reach Mach 7 and was attempted on June 2, 2001. The HXLV lost control shortly after launch, 
resulting in the loss of both the HXLV and HXRV (ref. 2). During Flight 2, flown to Mach 6.8 
on March 27, 2004, the HXRV successfully demonstrated the in-flight operation of the scramjet 
(ref. 3). All of the goals for that mission were achieved, including positive acceleration of the 
vehicle by the scramjet. The third and final flight was flown to Mach 9.6 on November 16, 2004 
(ref. 4). During both successful missions, the HXRV was in a controlled autonomous flight from 
the point of separation to splashdown in the Pacific Ocean. 

The Hyper-X Research Vehicle Mission Description 

The HXRV mission timelines were similar for Flight 2 and Flight 3. The principal difference 
between the two flights is the higher Mach number attained during Flight 3, which resulted in a 
correspondingly longer descent trajectory. Figure 2 shows an overview of the HXRV mission. The 
launch, scramjet experiment, and descent portions of the mission were conducted off the coast 
of southern California. For both missions, the HXLV propelled the X-43A stack to a separation 
altitude of approximately 100,000 ft and a dynamic pressure of approximately 1000 psf. The HXRV 
separated at Mach 6.9 during Flight 2 and at Mach 9.7 during Flight 3. Several seconds after 
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separation, the primary experiment, which consisted of the scramjet engine test, was conducted. 
The primary experiment lasted approximately 45 s for Flight 2 and 30 s for Flight 3. After the 
completion of the primary experiment, the cowl door was closed and a recovery maneuver was 
initiated to arrest the vehicle descent and dynamic pressure buildup. After the descent was arrested, 
an unpowered trajectory was flown to a splashdown into the Pacific Ocean. During the descent, a 
series of parameter identification (PID) maneuvers were performed at every integer Mach number 
down to Mach 2 (refs. 5, 6). These maneuvers began at Mach 5 for Flight 2 and at Mach 8 for 
Flight 3. 
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Figure 2. Hyper-X flight trajectory. 


MONTE CARLO ANALYSIS OVERVIEW 

The DFRC 6-DOF nonlinear simulation of the HXRV was used in conjunction with a series 
of Matlab® (The Math Works, Natick, Massachusetts) scripts to conduct the Monte Carlo analysis. 
The primary objectives of the Monte Carlo uncertainty analysis were to assess vehicle stability in 
the presence of modeling and measurement uncertainties, the flight conditions during the engine 
test, the unpowered descent trajectory, and the controlled flight impact area. Statistical information 
from the Monte Carlo analysis was examined, but was not the primary product. The project was most 
concerned with assessing the vehicle performance and the expected flight envelope in the presence 
of the expected model uncertainties. This performance assessment turned out to be largely Boolean 
in nature. The vehicle performance was either acceptable or not. Software performance, including 
the flight control laws, guidance, and scramjet controller (ref. 7) was a part of this assessment. The 
Monte Carlo objectives were achieved by conducting many thousands of simulation runs, with 
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each run having a random combination of modeling uncertainties. To gain a greater understanding 
of the relationships between the dispersions and their effects on the trajectory, data were correlated 
to determine combinations of factors that most influenced the performance of the aircraft. A 
secondary objective consisted of using the uncertainty analysis data to find a limited set of runs 
that could be used for thorough validation testing of the flight hardware and software systems. 

The Hyper-X Research Vehicle Monte Carlo Analysis Tool (MCAT) 

A Monte Carlo analysis tool (MCAT) was developed to stress the vehicle system with 
predicted uncertainties. The HXRV MCAT tests the vehicle performance in the presence of random 
and unique sets of input parameter uncertainties. The MCAT consists of a set of Matlab® M-files 
that are used as controlling scripts to drive the nonlinear simulation with various dispersions, 
analyze the results, save the data, and plot the data for further analysis. Matlab® M-files are used 
to automate the process of generating the thousands of simulation scripts. 

An execution of the MCAT consists of a large number of simulation runs with the first run 
being nominal and the rest containing different combinations of uncertainties. 

Manipulation and storage of the Monte Carlo results presented numerous challenges. Several 
gigabytes of data were produced by each set of simulation runs. The analysis tools run in the 
Matlab® workspace were incapable of handling a complete Monte Carlo dataset. The large amount 
of data resulted in the partitioning of the Monte Carlo data into multiple .mat files. A separate .mat 
file was generated for each simulation output variable and this file contained data from all of the 
simulation runs. To minimize the amount of data stored from a Monte Carlo analysis, trade-offs 
were made between the rate at which the data were saved and the number of signals saved for 
short-term and long-term storage. The parameters of greatest interest were saved at 100 Hz during 
the engine experiment phase. During other phases of the simulated missions, data were saved 
at 10 Hz. Relatively slowly-varying parameters, such as altitude, Mach, latitude, and longitude 
were stored at 1 Hz. Table 1 shows the Monte Carlo parameters saved at each rate. Parameters of 
greatest interest were retained indefinitely, whereas less important data was deleted after a week to 
allow storage space for the next set of Monte Carlo results. 

Figure 3 is a flowchart graphically representing the stochastic MCAT process. The inputs of 
each simulation mn, including the desired uncertainties, are defined in an input file that is called a 
simulation input script. The MCAT generates the simulation input scripts and details the number of 
simulation runs. The MCAT then uses the simulation input scripts to run the simulation in a batch 
mode a number of times specified by the MCAT and saves a subset of the simulation parameters 
that are recorded for later analysis. The simulation input scripts were deleted after each Monte 
Carlo run to save space. The simulation output data of each of these saved parameters are collected 
from each simulation run into separate matrices (one for each parameter). The data matrices are 
then stored for future analysis and plotting. 
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Figure 3. Monte Carlo analysis tool overview. 


For the MCAT, every input parameter that is varied has an uncertainty expressed in terms 
of a probability density function. The probability density functions for each input parameter are 
expressed in terms of its distribution type (uniform or Gaussian), a mean, and a standard deviation. 
The vast majority of the simulation input parameters are varied independently. Several parameters, 
such as the moments of inertia, were checked for consistency. If a parameter was found to be 
inconsistent with the other parameters with which it is correlated, the value was thrown out and 
recast until a consistent value was found. The input parameter values are based upon their particular 
probability density functions. They are selected within the ± 3 sigma bounds for each run with 
either a uniform or Gaussian distribution (the distributions are specified by the user in the MCAT). 
If a distribution is Gaussian, the 3 sigma bounds encompass 99.87 percent of the possible values 
of the uncertain parameter. The uncertainties associated with a parameter were provided by each 
discipline responsible for that portion of the model (i.e. aerodynamics, propulsion, flight control, 
etc.) (ref. 7). 




























































As a part of each Monte Carlo analysis, summary statistics were extracted from the dataset 
and placed in a Monte Carlo report. A temporary 100 Hz dataset was created to extract the statistical 
information from the highest simulation data rate available. The temporary dataset was deleted 
after the statistics had been gathered. A number of statistics, including the maximum and minimum 
values, were recorded for certain parameters during particular flight phases. In addition, the run 
number associated with each maxima or minima was recorded. Table 2 shows a sample page from 
the Monte Carlo report. This report was examined and compared to previous reports to determine if 
any unusual cases or events had been generated in the latest set of Monte Carlo runs. Effective use 
of these summary statistics helped to reduce the long-term data storage requirements and identify 
runs of interest. The simulation scripts for each Monte Carlo analysis could be easily regenerated. 
This allowed for a simulation mn to be repeated and more data acquired if, at a later date, a review 
of the summary or lower rate data indicated an item of interest. 

In preparation for the second Mach 7 flight attempt, 1600 Monte Carlo simulation runs were 
conducted with the latest models and software as part of the final analysis. Many thousands of 
Monte Carlo simulation runs had been conducted prior to the generation of the final dataset as 
the analysis had been run continuously as the simulation models had been refined and the vehicle 
software matured. The number of simulation runs used for the final analysis was determined by 
examining the statistics of notable parameters, such as the maximum value for filtered normal 
acceleration, seen in figure 4. The asymptotic trend of the parameters showed that there would 
be little expected change in the output of the Monte Carlo analysis when performing more than 
1600 runs. 


5 


4 

Mean and 
extremes, 3 
G 

2 

1 

0.8 

0.6 

Standard 
deviation, 0.4 
G 

0.2 

0 

0 200 400 600 800 1000 1200 1400 1600 

Number of runs 

070267 




Figure 4. Statistics for maximum filtered normal acceleration. 


8 
















SIMULATION DESCRIPTION 


The Hyper-X program developed a high-fidelity 6-DOF nonlinear simulation model of 
the HXRV. The HXRV simulation does not model the HXLV and was usually initialized to the 
separation conditions from the HXLV. A separate 6-DOF simulation was developed to model the 
HXLV, and is beyond the scope of this report. A simplified model of the separation dynamics 
is included in the HXRV simulation. A separate, more detailed simulation, known as SepSim, 
was developed to fully model the HXRV separation from the HXLV and is beyond the scope 
of this report. The HXRV simulation was based on the standard DFRC simulation architecture 
(ref. 8) and contained detailed models that included the vehicle aerodynamics, mass properties, 
atmospheric, and wind models. Vehicle subsystems modeled in the simulation included the flight 
control system, actuator, engine and fuel systems, the inertial navigation system (INS), and the flush 
airdata sensing (FADS) system. The flight management unit (FMU) software models consisted of 
guidance, navigation, and flight control systems, in addition to the propulsion system control. 
Since most of the components in the research vehicle were single-string (not fault redundant), no 
failure modes were modeled in the simulation. The majority of the simulation core source code 
and models are written in Fortran 77 with several of the models consisting of autocode generated 
from Simulink® (The Math Works, Natick, Massachusetts) models. The simulation graphical user 
interface is written in C programming language. 

The simulation can be run in either a real-time mode or a non-real-time batch mode. The batch 
mode operation is faster than real-time and was used for most engineering analysis including Monte 
Carlo analysis. The real-time mode was used only when the simulation was connected to ground 
test hardware. The base simulation frame rate is 200 Hz. Vehicle models and simulation routines 
are called at the appropriate multiples or submultiples of the base frame rate. Numerical integration 
of the simulation state information is performed using a second-order fixed-step Runge-Kutta 
algorithm optimized for fast execution. Oblate-Earth, 6-DOF equations of motion are employed. 
The simulation employs a command line interface along with a scripting capability. The scripting 
capability allowed for the simulation execution of script files containing multiple commands that 
had been strung together. The MCAT made extensive use of the simulation scripting capability 
to generate a unique simulation setup script for each Monte Carlo simulation run. Nearly every 
variable in the simulation is capable of being recorded in an external data file. The selection of 
output signals was done through the script files. 

MODEL UNCERTAINTIES 

The Monte Carlo analysis consisted of applying uncertainties to the simulation vehicle 
models. A total of 286 model uncertainties were used as a part of the final analysis for the second 
and third flights of the HXRV. The uncertainty distributions for each model were estimated using 
measurement error buildup, test data spread, or engineering judgment. The model uncertainties 
were generally developed by the engineering discipline team responsible for producing the 
model. These model uncertainties could be grouped in the following categories: separation 
condition, mass properties, control surface actuator performance, HXRV aerodynamics, separation 
aerodynamics, separation piston performance, FADS, atmosphere, winds, HXRV sensors, engine 


9 


performance, and propulsion system performance. Uniform uncertainty distributions, described by 
their minimum and maximum values, were generally used as they produced the most conservative 
results. Uniform distributions were used for separation conditions, engine performance, winds, and 
inertial navigation system (INS) errors. To describe each Gaussian distribution, four numbers were 
used: mean, standard deviation (1-sigma), minimum, and maximum. The minimum and maximum 
numbers were used to limit the dispersions to better match their intended distributions and to 
maintain a 3-sigma limit on the applied uncertainties. The Matlab® rand or randn function was used 
to generate a random number with either a uniform or Gaussian distribution for each uncertainty. 
These random numbers were then linearly mapped to the model uncertainty distribution to generate 
a model uncertainty parameter for input to the simulation. 

An overview of the Monte Carlo uncertainties is presented in Table 3. A detailed discussion of 
each area follows, with details presented where the data is not proprietary, sensitive, or otherwise 
restricted in distribution. 

The HXRV simulation was initialized at the point of separation from the HXLV for the Monte 
Carlo analysis. The flight conditions used to initialize the HXRV simulation were derived from the 
project requirements document and the HXLV Monte Carlo simulation results. The HXLV Monte 
Carlo results generally fell within specified bounds. When a specific project requirement existed, 
that bound was used to determine the mean and ± 3 sigma bounds for a particular condition since 
the requirement was more conservative than the HXLV Monte Carlo results. When no preexisting 
project requirement existed, the output of the HXLV Monte Carlo simulation was used. Table 4 
details the separation conditions and uncertainties used to initialize the HXRV simulation for the 
Monte Carlo analysis. 

The mass properties uncertainty model was developed by the HXRV structures team and is 
largely based on the accuracy of the test data. The vehicle weight, center of gravity, and principle 
moments of inertia were measured on the ground and adjusted for the final flight configuration. 
Measurement uncertainties and engineering judgment were used to determine the appropriate 
uncertainties. The mass property uncertainty values are presented in table 4. 

The HXRV control surface actuator models are discussed in detail in reference 9. A high- 
fidelity, proprietary model of the HXRV control surface actuators was provided by the manufacturer, 
Moog (New York) and updated by the HXRV guidance, navigation, and control (GNC) team, which 
also developed an uncertainty model for the actuator model. This high-fidelity model was used 
for the HXRV Monte Carlo analysis for Flights 2 and 3. The HXRV GNC team also developed a 
simplified second-order model that approximated the actuator performance. This simplified model 
was used for Monte Carlo analysis prior to the first flight attempt and these model uncertainties 
are presented in table 4. In general, the simplified actuator uncertainty model is more conservative 
than the high-fidelity model. 

The HXRV vehicle aerodynamics model was developed from wind tunnel testing and 
computational fluid dynamics (CFD) (ref. 10). The aerodynamic uncertainty model was developed 
based on the wind tunnel testing and CFD results (ref. 11). In addition, work done in developing the 
X-33 uncertainty model (ref. 12) was leveraged in the uncertainty development. The aerodynamic 
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coefficient uncertainty is a function of Mach number and engine operation. The aerodynamic 
uncertainties for cowl-closed flight decrease as the Mach number decreases. The aerodynamic 
uncertainties are increased for cowl-open flight and increased further while the scramjet is operating. 
The uncertainties in the principle aerodynamic coefficients during cowl-closed flight are presented 
as percents in table 4. Aerodynamic coefficient biases were used, but are not presented as they 
were not formulated in terms of a percent. The magnitude of the aerodynamic coefficient biases is 
generally small relative to the coefficients. 

A simplified separation aerodynamics model was developed which applied delta force 
and moment increments to the free-flight HXRV aerodynamic coefficients to approximate the 
aerodynamic interactions expected to occur during the separation event. This model is known as 
the Sep Delta model. Uncertainties on these delta increments were developed and incorporated 
into the Monte Carlo analysis. The Sep Delta model was derived from SepSim analysis. A detailed 
presentation of the Sep Delta model along with its uncertainties is beyond the scope of this report. 
The principle purpose of the Sep Delta model in the HXRV Monte Carlo simulation was to replicate 
the spread in flight conditions the vehicle could be expected to encounter after the separation event. 
A detailed discussion of the SepSim analysis and results is contained in reference 13. 

Two pistons were used to push the HXRV away from the HXLV during the separation event. 
A piston force model was developed, along with uncertainties, primarily from test data. A piston 
performance model was implemented in the HXRV simulation, in addition to SepSim. A detailed 
discussion of the piston performance model and uncertainties is beyond the scope of this report. 

A flush airdata sensing (FADS) system was employed on the HXRV (ref. 14). The FADS 
system provided an angle of attack estimate that was blended with the inertial measurement 
following the engine test during the Mach 7 flight. The FADS system was only used real-time 
during the unpowered descent portion of the Mach 7 flight. The FADS performance uncertainties 
were examined by introducing uncertainties in the simulated pressure port measurements. There 
are three components to the FADS pressure port uncertainties. White noise was added to each 
frame as a uniformly-distributed number continuously varying from -1 to 1. There was also an 
uncertainty in the pressure coefficient at each port, which was a function of Mach and alpha. 
Finally, a PPT (pressure port transducer) measurement uncertainty was added as a normally- 
distributed number between -1 and 1. For the Monte Carlo analysis, uncertainties were applied 
only to the pressure ports used as part of the angle of attack estimation algorithm. Table 4 details 
the FADS uncertainties. 

The HXRV simulation contains a separate atmosphere model for each month. This model 
was developed by the DFRC meteorological group and tailored for the expected HXRV engine 
test location. The model is derived from Global Reference Atmospheric Model (GRAM-95) 
data (ref. 15). The atmosphere model used by the HXRV onboard software is constant during a 
Monte Carlo run and is set to the expected flight month. Atmospheric uncertainties are modeled 
by choosing a uniform random month for the atmosphere model in the nonlinear simulation. By 
randomly varying the simulation atmosphere model between the 12 available models, adequate 
uncertainty was applied to the HXRV onboard software relative to the simulated values. 
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North-South, East-West, and vertical wind components are independently varied every 
10,000 ft. North-South and East-West wind data was provided by the DFRC meteorology group. 
Vertical wind mean data is from GRAM-95, and also varies by altitude and month. Because of 
sparse statistical data for vertical winds over the ocean, a worst-case 3 sigma value of 10 ft/s was 
used. Table 4 provides sample wind uncertainties. 

A high-fidelity sensor model, including uncertainties, was developed to model the performance 
of the INS contained in the H-764 Honeywell flight control computer (Honeywell, Inc., Minneapolis, 
Minnesota). The INS package is capable of utilizing global positioning system (GPS) updates to 
aid the position and attitude solution, but the HXRV used the pure inertial solution throughout its 
mission, so GPS performance was not included in the sensor model. Information used to develop 
the high-fidelity sensor model is proprietary and therefore, the model details are not presented 
here. 


The HXRV simulation contains an engine database that approximates the steady-state scramjet 
engine thrust and moments, hereafter referred to as performance (ref. 16). The simulation engine 
database consists of maximum and minimum values as a function of flight condition and fuel-to-air 
ratio. Engine performance is varied for each run by setting an engine performance parameter that 
selects a point between the minimum and maximum values in the database. Uniformly distributed 
values between 0 and 1 were used to vary the engine performance. The engine performance 
parameter was set at the beginning of a Monte Carlo simulation run and not varied thereafter. A 
value of 0.5 corresponds to the midpoint of the maximum and minimum values and is considered 
the nominal engine performance. 

Uncertainties in the HXRV scramjet fuel system performance were modeled in the Monte 
Carlo analysis. Uncertainty in the fuel flow, igniter flow, and airflow were also modeled. 

Most of the applied uncertainties were independent of each other; however, some 
reasonableness checks were performed on correlated parameters. These included angle of attack, 
pitch attitude, and flight path angle, as well as altitude, velocity, and dynamic pressure. Moments 
of inertia were also evaluated to ensure that the sum of two principle axes would be greater than 
the third axis. 

In addition, the HXRV simulation random number seed was set by the Monte Carlo tool. The 
random number seed is used to initialize the random functions in the HXRV simulation that apply 
noise to the sensors and FADS pressure measurements. This explicit setting of the random number 
seed allowed for repeatability of a single Monte Carlo run. 

MISSION AND TRAJECTORY SENSITIVITY ANALYSIS 

Key relationships between uncertainties and performance were determined using a sensitivity 
study. In general, the project was most interested in determining the uncertainty combinations that 
would produce the maximum or minimum of a particular metric. For example, one metric was the 
vehicle positive acceleration caused by the scramjet engine operation. As a result, the maximum 
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and minimum positive acceleration values were of interest along with any key uncertainties 
influencing the value. 

Performance metrics were unique to the various mission segments. The following were some 
of the metrics that were tracked during the separation event and parameter identification maneuvers: 
surface command rates and positions, vehicle attitude rates and angles, and normal acceleration. 
The most important performance metrics were the engine-on flight conditions: primarily angle 
of attack, Mach number, and dynamic pressure. These metrics were used to provide the scramjet 
team with the expected flight conditions for their own detailed off-nominal analysis of the scramjet 
performance. The robustness of the guidance algorithm to off-nominal engagements was assessed 
by examining the impact of off-nominal headings when the algorithm was engaged. Over the rest 
of the descent to the ocean, surface commands and positions, altitude error, vehicle attitude, and 
flight condition were tracked. The vehicle orientation, velocity, splashdown point, aim point error, 
and time aloft were then determined at the end of each run. This final information was used to 
evaluate the performance of the guidance routine and to develop a splash point boundary for range 
clearance. 

Maximums and minimums of these parameters of interest were drawn on scatter plots versus 
each model uncertainty input for a set of Monte Carlo runs to determine the uncertainty parameters 
to which the vehicle performance was most sensitive. An obvious example is seen in figure 5, 
which shows how increasing the drag (axial force aerodynamic coefficient) decreases the distance 
traveled. 
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mi 


500 



250 - 1 - 1 - 1 - 1 - 1 - 
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Figure 5. Hyper-X Research Vehicle down range versus drag dispersion. 
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The HXRV was found to be sensitive to a model uncertainty if it caused the vehicle to reach 
a maximum or minimum limit of some parameter of interest, such as angle of attack, normal 
acceleration, or control surface deflection. In addition, if the vehicle exhibited an undesirable 
behavior, such as a pitch or roll oscillation, the cause(s) were noted as parameters of interest. 
Before the first flight attempt, a ranked list was formed of the model uncertainties to which the 
HXRV performance metrics and trajectory were most sensitive. 

1. Separation angle-of-attack 

2. Separation sideslip angle 

3. Separation pitch attitude 

4. Separation heading 

5. Separation altitude 

6. Separation Mach 

7. Lateral center of gravity location 

8. Actuator damping ratio 

9. Axial force aerodynamic coefficient (approximately equivalent to Drag) 

10. Normal force aerodynamic coefficient (approximately equivalent to Lift) 

11. Pitching moment aerodynamic coefficient 

12. Error in pitch attitude caused by mounting and measurement 


Stress testing was conducted with combinations of these parameters being set at or beyond 
their minimum or maximum expected values. This testing resulted in a handful of cases that 
stressed the vehicle performance at varying points throughout the mission. These limited numbers 
of cases could be mn quicker than a full Monte Carlo analysis to assess vehicle performance, 
while still stressing the vehicle software and performance. These stress tests were used to quickly 
examine design changes to the flight control laws and guidance algorithm. In addition, the stress 
test analysis partially accounted for the unknown unknowns present in every flight research 
program by examining the vehicle performance with what are believed to be unrealistic model 
uncertainties. These tests provided an assessment of the vehicle robustness in excess of that which 
is expected to be required. 


VALIDATION TESTING 

Two types of flight hardware and software validation tests were performed: hardware-in- 
the-loop (HIL) and aircraft-in-the-loop (AIL) tests (ref. 17). The purpose of HIL testing was to 
validate the hardware and software performance of the flight control computer, which is a portion 
of the flight management unit (FMU). For HIF testing, the FMU was connected to the simulation. 
Simulated inertial sensor data was provided to the FMU by the simulation. The flight software 
running on the FMU generated control surface commands and propulsion control valve commands 


that were in turn fed to the simulation. The number of test cases run in the HIL environment 
was limited because of the test complexity and time involved. The AIL tests were similar to the 
HIL tests but also included the vehicle flight actuator controller and actuators as part of the test 
configuration. The AIL testing placed flight hardware at risk, so an even more limited number of 
test cases were chosen for this type of testing. 

Taking these factors into consideration, a small set of runs was chosen for validation testing 
from the thousands of Monte Carlo runs. Approximately 20 runs were performed using the HIL 
setup, and less than 5 runs were conducted as a part of AIL testing. The runs chosen for HIL and 
AIL testing were chosen from a Monte Carlo dataset by finding the minimum number of runs 
that bounded some desired aspect of the mission. For example, figure 6 shows the 5 cases found 
which bounded 90 percent of the altitude profiles from 1600 Monte Carlo cases. In general, the 
HIL cases were chosen to bound the flight conditions seen in the Monte Carlo analysis and the 
AIL cases were chosen to stress some aspect of the control surface actuators such as response to a 
large step change, maximum/minimum deflection, or high levels of surface motion. An automated 
tool was developed which selected the minimum number of cases that met the desired criteria. 
Table 5 shows the criteria the tool maximized or minimized during a given mission segment. The 
tool worked its way through the Monte Carlo dataset and selected a run for HIL testing if it was 
90 percent of the way to a specified bound. Cases exceeding the 90 percent threshold for multiple 
specifications were ranked higher and the highest ranked cases were chosen for HIL analysis. 


Altitude 



Increasing time since separation 
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Figure 6. Validation test cases compared to full Monte Carlo results. 
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FLIGHT DATA COMPARISONS TO PREFLIGHT MONTE CARLO RESULTS 


A summary of the HXRV in-flight performance comparison with preflight Monte Carlo 
results is discussed for the Mach 7 and Mach 10 flights, with an emphasis placed on highlighting 
modeling shortfalls. Reference 18 provides a detailed discussion of the HXRV performance and 
comparisons with Monte Carlo data during the Mach 7 flight and reference 19 provides the same 
for the Mach 10 flight. 

Flight results are presented for the time period after the separation event. The HXRV simulation 
only roughly approximates the separation event and this time period is fully modeled in detail by 
the HXRV separation simulation. 


The Mach 7 Mission 

For the Mach 7 flight, the in-flight performance, when compared to Monte Carlo predictions, 
indicates that the analysis was fairly accurate at predicting actual flight performance and bounds. 
Notable differences between the Monte Carlo predictions and the in-flight performance are 
discussed. Figure 7 displays the inertial angle of attack during the scramjet test. In figures 7-15, 
flight data is compared to a nominal simulation run started at the in-flight separation point. Monte 
Carlo data is represented on figure 7 by the gray bands that show the bounds for 66 percent of the 
runs, 95 percent of the runs, and the minimum and maximum values. 
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Figure 7. In-flight inertial angle of attack during scramjet test conditions compared to Monte Carlo 
predictions for the Mach 7 mission. 
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The cowl-open transient prior to the engine test resulted in the vehicle nose pitching up, 
whereas the nominal and Monte Carlo preflight analysis predicted a nose-down transient as shown 
in figure 7. The inaccurate prediction of the direction of the cowl-open transient was likely caused 
by an incorrect understanding of the unsteady aerodynamics caused by the cowl door opening. 

The inertial angle-of-attack response during the engine test is compared with the preflight 
Monte Carlo results in figure 7. During the engine test, the HXRV maintained the engine test 
conditions well within the requirements of ± 0.5° angle of attack. The inertial angle of attack fell 
outside the Monte Carlo bounds for a few seconds during the engine test as the igniter was removed 
from the fueling profile. This event is an example of the difficulties in modeling the complexities 
of the scramjet operation, as it was not covered by the Monte Carlo preflight predictions. The 
HXRV simulation contains a simplified engine model that is essentially a table lookup based on 
flight condition and therefore, does not fully model the dynamics of scramjet operation. 

Following the scramjet experiment, the HXRV increased the angle of attack to arrest the 
buildup in dynamic pressure and heating. During the recovery maneuver, the HXRV experienced 
small amplitude angle of attack oscillations at a frequency of approximately 0.65 Hz as shown in 
figure 8. The preflight Monte Carlo predictions did not include any cases with sustained oscillations 
during the recovery maneuver. The cause of the oscillations is not fully understood. After the 
recovery maneuver, the vehicle transitioned to the descent mode. 



—— Flight data 
—- Simulation 


Min/max 
95 percentile 
66 percentile 


070271 

Figure 8. Inertial angle of attack compared to Monte Carlo predictions for the Mach 7 mission. 
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The vehicle maintained controlled flight and performed well during the descent. Figure 8 
also shows that the inertial angle of attack was within the preflight predictions during the descent 
portion of the mission. The large changes in angle of attack during the descent are caused by 
preprogrammed PID maneuvers, which were preformed at integer Mach numbers. 

Figure 9 shows that the flight altitude profile was close to the nominal simulation results and 
was well within the Monte Carlo predictions. The bank angle was within Monte Carlo predictions 
throughout the descent except for a deviation late in the descent, seen in figure 10. At that point, the 
HXRV banked approximately 30° left wing down. The cause of this upset has not been conclusively 
determined, but the upset might have been caused by winds or unmodeled transonic aerodynamics. 
The impact location was well within the preflight prediction as shown in figure 11. Figure 12 
shows the HXRV ground track following separation from the HXLV overlaid on the preflight 
Monte Carlo predictions. The California coastline is included in figure 12 for approximate scale 
and orientation. 
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Figure 9. Inertial altitude compared to Monte Carlo predictions for the Mach 7 mission. 
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Figure 10. Bank angle compared to Monte Carlo predictions for the Mach 7 mission. 
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Figure 11. Mach 7 mission splashpoint compared to preflight Monte Carlo predictions. 
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Figure 12. Mach 7 mission ground track compared to preflight Monte Carlo predictions. 


The Mach 10 Mission 

The HXRV generally flew within the Monte Carlo bounds for the Mach 10 flight, with several 
notable exceptions (ref. 19). Figure 13 shows the in-flight angle-of-attack performance during the 
engine test. A slight pitch up, instead of a pitch down (as predicted by the nominal simulation), 
occurred when the cowl door opened. This pitch up occurred in both the Mach 7 and Mach 10 
missions. Before the Mach 10 flight, the engine database uncertainties were increased to account 
for the discrepancy seen in the Mach 7 flight. As a result, the in-flight angle-of-attack performance 
when the cowl door opened fell within the Monte Carlo predictions for the Mach 10 flight. 

The expected longitudinal upsets, when the engine was turned on and off, were within the 
Monte Carlo predictions. During the remainder of the engine test, the FIXRV maintained the required 
angle of attack. As seen in the Mach 7 flight, the inertial angle of attack fell outside of the Monte 
Carlo bounds during the engine test as the igniter was removed from the fueling profile. Figure 14 
shows that during the engine test, the AMW deflections were near the Monte Carlo bounds, which 
is believed to result from an inaccurate prediction of the nominal Cmo (pitching moment at zero 
degree angle of attack) by the aerodynamic and engine databases. This discrepancy indicates that 
a larger Cmo uncertainty might have been appropriate. 
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Figure 13. In-flight inertial angle of attack during scramjet test conditions compared to Monte 
Carlo predictions for the Mach 10 mission. 
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Figure 14. In-flight elevator position compared to Monte Carlo predictions for the Mach 10 
mission. 
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Following the engine test, the FIXRV successfully closed the engine cowl door and pulled to 
the desired angle of attack as part of the recovery maneuver. After approximately 10 s, the FIXRV 
began to slowly rock approximately ± 10° about the commanded bank angle, and the rocking grew 
to as much as ± 20° about the command when guidance started maneuvering the vehicle in the 
direction of the aim point. During the same time period, the internal engine pressures indicated 
that air had started to flow through the engine, resulting in a pressure increase on the aft underbody 
of the vehicle, which pushed the nose down. The FIXRV responded by commanding the elevators 
to move to a more negative deflection to bring the nose back up. Figure 15, which shows the left 
AMW deflection from flight compared to the Monte Carlo data, illustrates the primary effect of air 
flowing through the engine. The AMW deflection was outside of the Monte Carlo bounds at high 
Mach numbers. Another event of note during this time period was a step change in the elevator 
command of several degrees to counter a longitudinal disturbance. This disturbance is believed 
to have been caused by a venting in the fuel system and the spontaneous ignition of igniter in 
the engine during this process. The internal engine pressures dropped shortly before the Mach 8 
PID maneuver, and the HXRV flew generally as predicted from Mach 8 and below. The AMW 
deflections continued to be outside of the Monte Carlo bounds until approximately Mach 5 when 
they started to track the nominal preflight predictions. Unexpected airflow through the engine and 
the unexpected spontaneous fuel ignition during the venting process are examples of phenomena 
that were not anticipated and thus not modeled as part of the Monte Carlo analysis. It is not known 
if the AMW deflections at the higher Mach numbers were outside the Monte Carlo bounds solely 
because of air flowing through the engine or if a mismatch in the predicted aerodynamics is also a 
cause. The FIXRV flew as predicted during the rest of the descent after about Mach 5. 
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Figure 15. Hyper-X Research Vehicle Mach 10 mission elevator position compared to Monte Carlo 
data. 











OBSERVATIONS AND LESSONS LEARNED 


It is fairly obvious that a vehicle simulation is only as good as the models of which it is 
composed. It is also obvious that a Monte Carlo analysis is only as good as the uncertainties used. 
In the case of the HXRV, some sets of uncertainties were known better than others. For example, 
the mass properties of the vehicle were measured directly and therefore were very well known 
with the uncertainty being largely driven by measurement accuracy. On the other hand, the HXRV 
was the first airframe-integrated scramjet-powered vehicle operated in flight. As a result, there was 
greater uncertainty in the aerodynamics, cowl motion transients, and engine operation because of 
the limited flight experience with this class of vehicle and the limited amount and applicability of 
wind tunnel testing. These models carried correspondingly larger uncertainties and, based on flight 
results, possibly should have carried larger uncertainties. 

The aerodynamic and engine performance uncertainties used in the HXRV simulation 
were primarily applied as a scalar to the vehicle forces and moments. The flight data indicated 
inaccurate predictions in the vehicle Cmo in several instances; therefore an improvement to the 
HXRV aerodynamic uncertainty model would have been the addition of larger biases to the total 
vehicle forces and moments. 

Monte Carlo analysis is very useful for predicting vehicle performance in the presence of 
known uncertainties. If, however, the uncertainties are believed to be conservative, it can provide 
a false sense of security and the perception that the vehicle will be robust to anything that could 
happen. Monte Carlo analysis also cannot account for unknown unknowns. The HXRV Mach 7 and 
Mach 10 missions experienced several events that were not anticipated and hence not modeled in the 
simulation. The unexpected pitch-up when the cowl door opened and the oscillations that occurred 
during the recovery maneuver during the Mach 7 mission are excellent examples of unanticipated 
events not modeled by the simulation and therefore not a part of the Monte Carlo analysis. The 
HXRV GNC team attempted to account for these unknown unknowns by maintaining conservative 
uncertainty bounds, primarily using uniform uncertainty distributions to produce a larger number 
of cases at the “corners,” and performing stress tests past the uncertainty bounds. Although the 
exact uncertainties and unknown unknowns were not captured within the Monte Carlo bounds, the 
system was robust enough to account for them without catastrophic consequences. 

The Monte Carlo analysis tools turned out to be fairly complicated. Keeping track of the 
hundreds of model uncertainties and the interactions between the Matlab®-based tools and the 
simulation turned out to be a straightforward, but involved process. The simulation models and 
uncertainties evolved greatly between the initial Monte Carlo capability developed for Flight 1 and 
that used for Flight 3. It is extremely important that any Monte Carlo analysis tool be straightforward 
to use and easily modifiable, especially when the tool is projected to be in use for many years. 

The Monte Carlo analysis tools were not under configuration control until late in the project. 
This lack of configuration control was not a serious problem since modifications to the tool were 
typically done by only one or two people working closely together, however, the HXRV simulation 
was under configuration control from the beginning. It became desirable to place the tool under 
configuration control later in the program since it was used for an important preflight analysis. The 
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majority of the uncertainty values were contained in the Matlab® scripts, which were a part of the 
Monte Carlo analysis tools, however, a number of the uncertainty values were contained in the HXRV 
simulation and a random number was input from the tools to set the uncertainty in the simulation. 
There would have been several advantages in placing the uncertainties in the HXRV simulation. 
First, all of the uncertainties would have been in the same tool and been under configuration 
control from the beginning. Second, the Matlab® scripts that generated each simulation run could 
have been much simpler if they had only to generate a set of random numbers for input to the 
simulation. An alternative approach would have been to incorporate all of the Monte Carlo run 
generation capabilities in the HXRV simulation. The main drawback to either of these approaches 
is that the simulation would have grown even more complex and potentially difficult to modify. 
The decision to place the Monte Carlo uncertainties and capabilities in the vehicle simulation is 
one that is specific to each project and depends on the complexity of both the simulation and the 
desired Monte Carlo analysis. 


CONCLUSION 

An overview of the Hyper-X research vehicle Monte Carlo analysis tools, model uncertainties, 
and analysis was presented. The methodology for selecting hardware validation test cases from the 
Monte Carlo data was described. The results of the preflight Monte Carlo analysis were compared 
with flight data, with an emphasis placed on highlighting modeling shortfalls. Modeling and Monte 
Carlo analysis lessons learned were presented. The preflight Monte Carlo analysis indicated that 
the X-43A control system was robust to the preflight uncertainties and provided the Hyper-X 
project an important indication that the vehicle would likely be successful in accomplishing the 
mission objectives. The X-43A accomplished all of the objectives during the two successful 
flights, however, the actual in-flight performance was not entirely captured within the Monte Carlo 
bounds. The X-43A flight control laws and guidance algorithms were robust enough as a result of 
the preflight Monte Carlo analysis, to ensure that the vehicle was able to handle the unexpected 
without undue consequences. 
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Table 1. Monte Carlo parameters saved at 1, 10, and 100 Hz. 


Description 

1 Hz 

10 Hz 

100 Hz 

Angle of attack 


X 

X 

Sensed angle of attack 


X 

X 

Angle of attack command 


X 

X 

Estimated angle of attack 


X 

X 

Axial acceleration 


X 

X 

Side acceleration 


X 

X 

Vertical acceleration 


X 

X 

Inertial angle of sideslip 


X 

X 

Angle of sideslip 


X 

X 

Coefficient of drag 


X 


Coefficient of lift 


X 


Cowl door command 



X 

Cowl door brake 



X 

Left wing command 


X 


Right rudder command 


X 


Right wing command 


X 


Aileron command 


X 

X 

Elevator command 


X 

X 

Rudder command 


X 

X 

Right rudder deflection 


X 

X 

Left wing deflection 


X 

X 

Right wing deflection 


X 

X 

Estimated engine fuel to air ratio 



X 

Actual engine fuel to air ratio 



X 

FADS angle of attack bias 


X 

X 

FADS dynamic pressure estimate 


X 

X 

FADS valid flag 


X 

X 

Time since separation 



X 

Flight path angle 


X 

X 

Altitude 

X 

X 

X 

Rudder hinge moment 


X 


Left wing hinge moment 


X 


Right wing hinge moment 


X 


Reference altitude 

X 

X 


Normal acceleration command 


X 


Geodetic latitude 

X 


X 
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Table 1. Concluded. 


Description 

1 Hz 

10 Hz 

100 Hz 

Sensed latitude 


X 


Sensed longitude 


X 


Mach number 

X 

X 

X 

Sensed Mach number 


X 


Pressures, valve settings, etc. 



X 

Axial acceleration 



X 

Average axial acceleration 



X 

Axial acceleration at CG 



X 

Side acceleration 



X 

Average side acceleration 



X 

Side acceleration at CG 



X 

Normal acceleration command 


X 


Normal acceleration 


X 

X 

Average normal acceleration 



X 

Normal acceleration at CG 



X 

Roll rate 


X 

X 

Bank angle 


X 

X 

Bank angle command 


X 

X 

Sensed bank angle 

X 

X 

X 

PID complete flag 




PID start flag 




Cowl door position 



X 

Heading angle 


X 

X 

Pitch rate 


X 

X 

Dynamic pressure 

X 

X 

X 

Sensed dynamic pressure 


X 

X 

Yaw rate 


X 

X 

Range to go 

X 

X 

X 

Pitch attitude 


X 

X 

Time 

X 

X 

X 

North velocity 


X 


East velocity 


X 


Down velocity 


X 


Latitude 

X 

X 

X 

Longitude 

X 

X 

X 
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Table 2. Monte Carlo statistics from the beginning of the Mach 8 PID to splashdown. 


Run # 

Maximum 

differential 

596 

- 

596 

666 

951 

547 

547 

629 

629 

785 

629 

951 

907 

629 

Minimum 

925 

- 

951 

515 

535 

620 

620 

239 

455 

401 

904 

629 

620 

239 

Maximum 

751 

- 

951 

oh 

581 

668 

668 

<N 

<N 

284 

251 

668 

668 

546 


Maximum 

differential 

0.44 

o 

0.44 

CO 

O 

o 

0.22 

4.17 

4.17 

i/O 

O 

OO 

7.76 

0.72 

9.32 

1.24 

1.79 

2.15 

Standard 

Minimum 

CO 

oo 

o 

1.77 

0.24 

ZVO 

OO 

OO 

© 

0.88 

1.54 


8L1 

CO 

0.87 

08 

OO 

O 

vq 

Average 

Minimum 

2.25 

80 

2.3 

-0.77 

-1.14 

-9.21 

-9.21 

-0.59 

-0.33 

-4.14 

OO 

© 

1 

uo 

OO 

-9.12 

0.24 

Minimum 

oo 

80 

-2.35 

-1.22 

-2.87 

-11.96 

-11.96 

-5.33 

-5.11 

-7.5 

-3.93 

-10.06 

-11.98 

-4.53 

Average 

All 

6.54 

80 

6.58 

-0.36 

O 

© 

0.21 

0.21 

5.72 

°q 

uo 

ZVO 

OO 

uo 

910 

90'0 

5.99 

Standard 

Maximum 

2.58 

o 

2.5 

o 

© 

ZVO 

OO 

O 

1.08 

3.76 

3.64 

0.94 

4.45 

160 

1.08 

3.88 

Average 

Maximum 

10.3 

80 

zvo\ 

-0.21 


8.89 

8.89 

16.64 

16.77 

4.48 

16.4 

5.98 

8.54 

16.84 

Maximum 

15.01 

80 

16.67 

-0.07 

3.63 

12.22 

12.22 

20 

20 


23.9 

8.78 

11.96 

21.3 

Units 

deg 

deg 

deg 

o 

deg 

deg 

deg 

deg 

deg 

deg 

deg 

deg 

deg 

deg 

Parameter 

Angle of 
attack 

Angle-of- 

attack 

command 

Estimated 
angle of 
attack 

Axial 

acceleration 

Angle of 
sideslip 

Left rudder 
command 

Right 

rudder 

command 

Left wing 
command 

Right wing 
command 

Aileron 

command 

Elevator 

command 

Rudder 

command 

Right 

rudder 

position 

Left wing 
position 
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Table 2. Continued. 
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Table 2. Concluded. 
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Table 3. Monte Carlo uncertainty models and summary of parameters varied. 


Uncertainty model 

Number of 
parameters 

Parameters varied 

Uncertainty 

distribution 

Separation condition 

12 

Angle of attack, angle of 
sideslip, altitude, longitude, 
latitude, Mach number, 
attitude angles, angular rates 

Uniform 

Mass properties 

8 

Center of gravity position, 
moments of inertia, weight 

Uniform 

Control surface 
actuator performance 

66 

Actuator performance 
parameters for the 2 all- 
moving wings and the 

2 rudders 

Uniform 

HXRV aerodynamics 

17 

Aerodynamic coefficients 
and biases 

Gaussian 

Separation 

aerodynamics 

9 

Scale factors on separation 
aerodynamics deltas 

Uniform 

Piston performance 
(separation) 

7 

Piston performance terms 

Uniform 

Flush airdata system 
(FADS) 

8 

Pressure port CP and 
measurement uncertainty 

Uniform 

Atmospheric 

conditions 

1 

Month 

Uniform 

Winds 

12 

Altitude (10-120) kft 

Uniform 

12 

Fiorizontal speed 

Uniform 

12 

Wind direction 

Uniform 

12 

Vertical wind velocity 

Uniform 

Sensors 

51 

Gyro error and bias terms 

Uniform 

42 

Accelerometer error and bias 

terms 

Uniform 

12 

INS initialization and 
misalignment errors 

Uniform 

Engine performance 

1 

Engine performance scalar 

Uniform 

Propulsion 

performance 

4 

Propulsion performance 
parameters 

Uniform 


Total number of 


parameters 


32 























Table 4. Monte Carlo model uncertainties. 
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Table 4. Continued. 
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Table 4. Continued. 
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Table 4. Continued. 
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Table 4. Concluded. 
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Table 5. Validation test selection criteria. 

[X denotes a parameter was selected for a maximum or minimum value] 


Parameter 

Separation 

Engine Test 

Recovery 

Maneuver 

Descent 

Splashpoint 

Angle of 
attack 

X 

X 

X 

X 


Angle of 
sideslip 

X 

X 


X 


Left wing 
deflection 

X 

X 


X 


Right wing 
deflection 

X 

X 


X 


Right rudder 
deflection 

X 

X 
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Flight path 
angle 



X 
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Altitude 


X 




Altitude error 
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Mach 
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Normal 

acceleration 
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Roll rate 
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Pitch rate 

X 

X 




Dynamic 

pressure 
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Yaw rate 
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X 




Phi 
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Theta 
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Time 
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